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We propose a new type of Monte Carlo approach in numerical studies of quantum sys- 
tems. Introducing a probability function which determines whether a state in the vector 
space survives or not, we can evaluate expectation values of powers of the Hamiltonian 
from a small portion of the full vector space. This method is free from the negative 
sign problem because it is not based on importance sampling techniques. In this paper 
we describe our method and, in order to examine how effective it is, present numerical 
results on the 4 x 4, 6 x 6 and 8x8 Heisenberg spin one-half model. The results indicate 
that we can perform useful evaluations with limited computer resources. An attempt 
to estimate the lowest energy eigenvalue is also stated. 
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1. Introduction 

A great many methods have been investigated so far to numerically calculate various quantities 
in quantum spin systems. They are classified into two categories, the exact diagonalization^^ and 
the Monte Carlo^'^-* approaches. Both of them provide us very useful ways to study quantum spin 
systems but, as was repeatedly reported, the former is difficult to apply to large size systems and 
the latter often suffers from the negative sign problem. In this paper we introduce a new method 
based on the power method^' in order to evaluate eigenvalues of these systems. Our purpose is 
to show that one can effectively investigate the ground state by means of a stochastic selection of 
states. 

The method we propose here is a kind of Monte Carlo approach, where stochastic variables 
play an important role. It, however, differs much from the conventional quantum Monte Carlo 
methods which employ random walks or importance samplings. '^^ In our method random variables 
are used to reduce the number of states of the vector space which is huge for most systems of large 
sizes. To this end we consider a new type of the probability function to which we refer as on-off 
probability function in this paper. Since it ^''switches off^ many states in the vector space we can 
calculate approximate expectation values of powers of the Hamiltonian from a small number of the 
"on" states. Repeating this process we can obtain averaged values which are very close to the exact 
values. Physical quantities such as the energy eigenvalue of the system would then be estimated 
from these expectation values with some additional assumptions. Throughout this paper we use the 
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two-dimensional Heisenberg spin one-half model as a concrete example in numerical calculations. 

In the next section we define the on-off probability function and explain its analytical properties. 
We then show, in a simple instance using the exact ground state for the 4x4 Heisenberg spin system, 
how we make use of this probability function when we calculate a quantity concerning to the vector 
space. Numerical results are presented in order to demonstrate that it nicely works. Section 3 is 
to describe our approach in a quantum system whose Hamiltonian is H. We define random choice 
matrices and present a way to stochastically calculate expectation values of (V' | \ with 
them. Two cases of | ■0) which draw our interest are then discussed. One is that the state is the 
exact eigenstate of ^f, | V') =| V'e), and another is an approximate of | -^e) which we denote as 
I V'a)- In section 4 we present our numerical results on the Heisenberg spin systems. Results on 
a 4 X 4 lattice, where it is easy to obtain | ^e), are given in detail so that we can make a close 
inspection of the method. Then results on the 6x6 and 8x8 lattices are presented together with 
the description of the state | ^a) we generated. We also show several assumptions we employed 
to estimate the energy eigenvalues for these lattice sizes. The values we present in this section for 
the 6x6 and 8x8 lattices are compatible with the ones reported in the preceding studies with 
different approaches. The final section is devoted to summary and comments. 

2. On-off Probability Function 

Let us introduce a probability function which is given by 



where a is a constant which is greater than or equal to 1. In this paper we refer it as the on-off 
probability function because the variable -q takes only two values, a (on) or (off), with the 
probability 1/a and 1 — 1/a, respectively. Clearly 



where (()) denotes the statistical average. In order to demonstrate the role of this probability 
function we consider, as a simple example, the sum of absolute values of all (non-zero) coefficients 
in a state | V') = S I O^i (q / 0) of the two-dimensional Heisenberg spin model, 



Pir,) = ^d{r,-a) + il-^)6ir,) 



(1) 




(3) 



(2) 



(4) 



n 
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and try to evaluate it by summing up less number of |ci|'s. For this purpose we calculate the 
following sum 

n 
i=l 

where {rj} denotes a set of random variables {r]i,r]2,--- ,??n} generated by the following on-off 
probability functions 

Piim) = -S{rii-ai) + {l~-)6{r)i), (7) 
tti being max{\,e/\ci\) for a given constant e. It is easy to see that 

{{S{r,})) =Y.\Cr\{{m)) = S. 
i=l 

The variance of S^j^y is given by 

= {{Sy) - {{Si.^r = E ieh\ - IqP) (8) 

|Ci|<6 

where {{r]ir]j)) = 6ij {{r]f)) + (1 - Sij) {{rn)) {{r]j)) = dijUi + (1 - Sij) is used. 

Let Ni = 1(0) if rii 7^ 0(= 0) and N = J2^i- Then we know by ((A^)) how many non-zero r)iS 
appear in (6). Since 

PN,{Ni) = -6{Ni-l) + {l--)6{Ni), 

f 00 



1 

{{Ni}) = / NiPNm)dNi = -, 

Jo O-i 



we obtain 



m = ((EA^.))= E 1+ E V- 

i=l \ci\>e \ci\<€ 

In numerical work we substitute the average over ngmpi randomly generated {i]}'s for the 
statistical average ((•••)), namely we measure 

where {r]}k denotes the A;-th set of the generated random variables {r/} = {rii,r)2,--- ,%}• The 
variance is then given by 

PS = (('S'{^}»smpl - {{S{r,}))smp\- (H) 

As the statistical error we employ 
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Note that, according to the Chebyshev inequahty in statistics, 



{{S{r,}))smp\- Er < < {{S{r,}))smpl + Er 

in 75 % confidence level provided that Ug can be replaced by pg.^^ 

The results for the exact ground state of the 4x4 Heisenberg spin model, which has 12870 non- 
zero coefficients, are shown in Table I. In this case it is straightforward to calculate S(= {{S^^y)) ), cjg 
and {{N)) from (5), (8) and (9) since we know all Cj of the state. In Table I we see that {{S^^j))smpi, 
Pg and ((A^))siupi calculated with ngmpi = 10^ are in good agreement with them. We also see that 
e = 0.2 is enough to estimate S in the precision of 0.1%. In other words we can estimate S in the 
precision of 0.1% from only 367 of 12870 non-zero coefficients of the state. 

3. Powers of Hamiltonian 

In this section we discuss on the expectation values of the L-th power of the Hamiltonian H. 
First we make a brief comment on the power method to evaluate the maximum eigenvalue of H. 
Suppose the eigenstates of H are | ^pi), H \ ipi) =\ tpi)Ei {i = 0, - ■ • , Ny — 1) where Ny is the size 
of the full vector space, and | |>| -E'i I for i > 0. In the power method one repeatedly operates 
H to & trial state | tp) , which is expressed as 

iVv-l 

I V') = 5Z I ^^)^- (13) 

so that one obtains for sufficiently large L 

H''\i^) = ^ \ i^i)hiE^ ~ I i^o)hoEl^ (14) 

1=0 

and 

(V' I i?^ I V) - hlE^. (15) 

Note that from the expectation values (V' | \ tl^) one can acquire information relative to the 
energy eigenvalue. For large size systems, however, it soon becomes impracticable to calculate 

I tp) or {ip I I ip) because one can not keep the huge number of coefficients for \ tp) in 
the full vector space. One therefore has to do with truncated vector spaces. 

Now we describe our formulation. Let us start with a state 

Nv 

1=1 

where {| i)} denotes an arbitrary basis of the full vector space. In our method we stochastically 
find a limited vector space where we can calculate the expectation values E{L) = {ip | | ip) for 
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large L. For this purpose we introduce an Ny x Ny diagonal matrix Mj^j. , which we call random 
choice matrix hereafter, 

/ r?i ••• ^ 



772 ■ • ■ 



\ ■■■ VNv J 



(17) 



with random variables rji determined according to (7). When Cj = we employ e/6 as aj in (7), 
where the parameter d{< e) is a positive constant independent of i. 

Using independent random choice matrices M^^(m)y = diag.{r]i^\r]2^\ ■ ■ ■ ,r]^^} {m = 
1, 2, • • • , L + 1) we define 

E^^}{L) = (V^|M{^(^+i)}i?M{^(^)}i?M{^(^_i)} • • • i?M{^(i)}|V^). (18) 

Here we use the coefEcient Cj of the initial wave function | il)) in the on-off probability function 
(7) for any m . In order to explain the essential point of the method in numerical studies, let 
nb(m) and na(m) the number of non-zero components of the state ^fM|^{m-i) j. • • • HM^^(i_)^\'il)) and 
M|^(m)j.-ff-M|^(m-i)} • • • -ffM|^(i)||V'), respectively. We first keep nb(l) non-zero components of | -0). 
Then we operate M|^(i)| to | V')- Since many of the random variables rjf^ (i = 1, 2, • ■ • , N\) are 
zero, the state M|^(i)i, | V') has much less non-zero components {i.e. na(l) ^ nb(l)). The number 
of non-zero components increases when we next operate H to M|^(i)j. | iji) (nb(2) > 7T,a(l)), but it 
drastically decreases after operating M|^(2)j to the state (na(2) <C nb(2)). Choosing the parameter 
e in the on-off probability functions we can repeat the operation up to L within the range of our 
computer facilities. 

Prom (18) we are led to, with hij = {i\H\j), 



i j k I 

i j k I 

= EE---EE^*^^^---^fc^'^' = ^^i^''i^) = ^(^) (19) 

i j k I 

because {{rji^^'nj^^ ' ' 'V^k'^Vi^'^^^)) — {{Vi^^)){{'nj^^)) ' ' ' {{v^k'^)){{Vi^~^^^)) — 1) which is guaranteed by 
the fact that r/^™^ and r/j™ are mutually independent for m ^ m' . Henceforth wc abbreviate 
{{'Hi^^)) ^iid (([^i™^]^)) {{rji)) and {{rjf)) because they do not depend on m. It should be noted that 
the variances of r/^™^ are the same for all m. 
The variance of E^jjj{L) is 
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j,j' k,k' 1,1' 

= X] X] " ' X] X] ^^^^ij ■ ■ ■ ■ ^i'^i'f ■ ■ ■ ^k'l'Cl' ■ {1 + Sii' ( )) - 1)} 

k,k' 1,1' 

X {1 + iiv])) - 1)} • • • {1 + 6kk'i {{vD) - 1)}{1 + Su'i {{ri)) - 1)} 
- E{Lf, (20) 

where we use {{rt'^vt^)) = (1 - 5n0((^.))((^^'» + h^' = 1 + ^n'( - 1). 

Let us add a few expressions which are useful to check our numerical values presented in the 
next section. They are obtained from (19) and (20) when the state | V') is an exact eigenstate 
I V-e) = E I V'i)/i {ft + 0) of namely when ii'lV'E) = \i>^)E or Y.j hijfj = Efi. 

{{EE{r,}{m = E^{L) = E\ (21) 
4{^}(1) = 2£;2^//(((r?|))-l) 

i 

+ E f'f'ihrjfi iiv!)) - 1)( iiv])) - 1), (22) 
al^^yi2) = 3£;^j;//(((r?f))-l) 

i 

+ 2E' Yl fiffihijfi {{vn - 1)( iiv])) - 1) 

i,k 

+ E fifk{h^Jfih,kf{ iivh) - 1)( ((^D) - 1)( ((^^)) - 1) , (23) 

where we added the suffix E to clearly show that the quantities are for an exact eigenstate. Note 
that we can analytically evaluate them using {{rjf)) = fli = e/ | /« |. 

In numerical study we first consider the case | tp) =| ip-E) on a small lattice, and then proceed 
to the case that the state is an approximate one denoted by | V'a)- In both cases we measure 
((£'{^}(L)))smpi and p^^y[L) from ngmpi samples in the same manner as (10) and (11) in the previous 
section, namely, 

(L)))3„,pi = 5^i?{,},(L), (24) 

pU(L) = (([^M(i^)]'))smpl- ((£^W(L)))smpl- (25) 
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The error of E^^y (L) is evaluated by 



EriL) = 2J-^^. (26) 

y '^smpl 



4. Numerical Results 



In this section we report our results for the two-dimensional Heisenberg spin model. The state 
on each site is represented by a 2 component of the spin. Let us first concentrate on the case of 
4x4 lattice. This lattice size is suitable for making a careful examination of the method because 
we can start with the exact ground eigenstate | ipE) and its eigenvalue E. Table II (III) shows the 
L = 1 {L = 2) results on ((£^E{77}(-^)))smpi obtained with 10^ samples for various values of e. Values 
of {L) calculated from (22) or (23) are also presented in the table. We see that cTe^^}. {L) and 
P^^^y{L) are in good agreement. We also see that the exact value E = —11.2285 lies in the range 
((EE{^}(l)))smpi ± Er{l) except for e = 0.05, while | ((£;E{r,}(2)))smpi - E'^ \ < Er{2) for all values 
of e. In Table IV we present ((£^E{77}(-^)))smpi up to L = 10 fixing e = 0.1 or e = 0.01. We observe 
that the error increases as L does, but its dependence on L is weak. The relative error for e = 0.1 
ranges from 0.11%(L = 1) to 0.48%(L = 10) while for e = 0.01 it is between 0.006%(L = 1) and 
0.015%(L = 10). Note that Er{L) decreases for smaller values of e, which however costs more 
computer resources. Figure 1 plots distributions of E-^^^y^{L) (L = 1,2,3,5,10), which shows 
the deviation from the exact 5 function. The abscissa is 50 times of the ratio of each E-^^^y^{L) 
{k = 1,2, ••• ,7T-smpi) to E-e{L). For small values of L we see that the distribution is symmetric 
about the exact value E'" and the peak is high. For large L, on the other hand, the peak becomes 
lower and the asymmetry grows; the data localize in the region where £'e{,,}^ (-^^) < E^ and the 
distribution has a long tail toward large value of £^E{T)}k (-^) • These features should be taken into 
account in the process of the error estimation. 

How about the number of states that appear in the calculation? We measure ((-/V))sinpi 
before and after we operate M^^(L)j to the state HM^^(L-i)y • • ■ HM^^(i)y\i/j-E), which we denote 
{{N^{L)))snip\ and ((A^f (L)))smpl, respectively Table V shows the results up to L = 10. For e = 0.1, 
we observe that {{N^{L)))smpi ^ ((^b'(^)))siiipi for each L and ((iV^(L)))smpi shrinks to almost 
a half of the whole vector space. This observation proves the merit of using the random choice 
matrices. When we set e = 0.01, however, the results are not so appealing because wc have to deal 
with almost all states of the vector space every time wc operate H. Wc also sec in the table that, 
as L grows, {{N^ (L))) s^^^pi converges to a constant which depends on the value of e. 

Let us show the results on the approximate state then. Here we employ |V'a) = Cn S|/^|>e 
for a small constant c, where Cat is the normalization factor. In other words wc abandon the 
coefficient fi of the exact eigenstate \iPe) if < c. For c = 5 x 10^'^ the number of non-zero 
component of IV^a) is 5382 and we obtain ('i/'Al-ff IV'a) = —10.7445, which should be compared to 
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(V'eI-H'IV'e) = -11-2285. Using this IV-a) we measure ((£;A{^}(i>)))smpi, ((Ar^(L)))si„pi and 

((Ar^(L)))sinpi up to L = 10 with ngmpi = 10^ and compare them to exactly calculated ('0A|-ff^|'0A)- 
Errors of Ea{^}(L)'s are evaluated by 2^ p\^^^{L) / Usrapi- Table VI and Fig. 2 show the results 
with 5 = 0.01. In Table VI we see a good agreement between ((£^A{?7}(-^)))smpi and (V'aI-H^^IV'a) 
for all L. We also see in Fig. 2 that for each L the shape of the distribution of £^a{^Jj^(L) is 
similar to that of £?E{7?}k(-^) ™ l^ig- 1 although the peak is slightly lower and the asymmetry is 
more outstanding. It should be noted that, as is expected, for larger L we obtain better evaluation 
for the true value Eq. In Table VI for example, we see Ea{1) = -10.7445 > - ^^^1(10) = 
- v^O.297136 X IQii = -11.1505 > = -11.2285. 

Finally we report our results on lattices larger than 4x4. Here we use 

E[^y{L) = I #M{^(^)}^M^^(^_i)} • • • ^M^^(i)} I i^). (27) 

instead of (18) in order to decrease the statistical errors. 

In order to apply our method to large systems we have to obtain \iIja) without knowing what 
the exact eigenstate is. Among several attempts we made, we took the following procedure 
which is based on the Suzuki- Trotter formula. 
( a ) Choose one trial state iV'trl)- 

( b ) Operate exp{-AHi) ■ exp{-AH2) to iV'trl), where Hi + H-z = H and A = 0.4. 
( c ) Abandon the small coefficients in exp{—AHi)-exp{—AH2)\iptrl) if necessary so that the number 
of non-zero components of the state is within the limit of our computers. 

For the 6x6 square lattice the number of non-zero components of the \iPa) thus obtained is O(IO^) 

with the minimum absolute value of the non-zero coefficients |c|7722n — 

1.14 X 10""^. 

Then we apply the on-off probability method setting e = 0.0057 and 5 = \c\min- Table VII 
shows values of {{EA{n}iL)))smph {{^h i^))) smpi and ((iV^(L)))smpi up to L = 5 with 5.5 x 10^ 
samples. The memory we need here is 356 MBytes including the hash table and the CPU time 
consumed to one sample is 1100 seconds by a Pentium III machine. 

Now we have obtained reHable estimations of Ea{L) = (V'a | \ iJja) (L = 1, 2, • • • ,5). How 
can we estimate physical quantities, Eq for instance, using them? We see that ^((£'A{,,}(5)))smpi = 
-23.68683 ± 0.00083 is lower than Ea{1) = -23.56093. But, because L is not large enough, even 
((£^A{r?}(5)))smpi provides only a poor upper bound for the true value Eq = —24.4394 which is 
obtained in ref.(5) by the exact diagonalization. One better way for the estimation would be to 
calculate the lowest eigenvalue E starting from the basis {|V'a), -f^lV'A), -f^^lV'A)}- It is easy to 
obtain the orthonormalized basis {iV'o)) iV'i)) IV'2)} from {\iPa),H\iPa),H'^\iPa)} and to calculate 
the Hamiltonian matrix elements which are the functions of Ea{L) {L < 5). It, however, turned 
out to be ineffective too; the obtained value is = —23.88, which is still too far from Eq. So 
we need additional assumptions to make a better estimate of the eigenvalue. Here, based on the 
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low energy property described by the linear spin wave theory or more sophisticated models/^ we 
empirically assume 



(V'A I I V'a) = qoE^ + qi f (xE)^ ■ x'^dx 

Jo 



= ^"'(^0 + TT^Vt) = ^(^' ^' '^O' ' (28) 
L + a + L 

where qq, qi and a are free parameters which should be determined together with E by the fit. Note 
that qo + qi/{a + 1) = 1 should hold due to the normalization of \tpA)- We look for the minimum 
of the difference D with -Lmax = 5, 

-^max r 1 1 TT' ( T \ \\ n2 



L=l 



^ _ ((-^A{)7}(-^)))smpl 

F(L,£;,go,9i,Q;) 



changing these four parameters and accept them if | go + <li/{oi + 1) — 1 |< 0.003 is fulfilled and D 
is less than the sum of the relative errors, which amounts 



E 

L=l 



Er{L] 



{{EA{r,}{L)))smpl 



8.5 X 10 



-9 



By this fit wc obtain —24.32 < E < —24.13, which indicates that the evaluation is useful. 

It is possible to advance to 8 x 8 lattice in the same manner, where the number of non-zero 
components of IV'a) is O(IO^) and |c|inin = 3.28 x 10~^. The results with n-smpl = 1-1 x 10'\ e = 0.01 



and 6 = |c|min are shown in Table VIII. Note that ^((-EA{,,}(5)))smpi = -41.03348±0. 03817 is lower 
than Ea{1) = -40.74998. By the fitting using (28) we obtain an estimation -43.5 <E< -41.1, 
which includes —43.103, the value also reported in ref.(5). It is large sample errors, which would 
be improved by increasing nsmpb that cause the wide range of the evaluated E here. 

5. Summary and Comments 

In this paper we have suggested an alternative Monte Carlo approach which, using the stochas- 
tic variables following the on-off probability function, does not rely on random walks nor importance 
samplings. Our purpose is to show the way to study large size systems without keeping the full 
vector space. We apply this method to the two-dimensional Heisenberg spin one-half model, which 
was chosen as a familiar example, and carry out numerical evaluations for lattices up to 8 x 8. The 
results show that we can obtain reliable numerical data on expectation values of the power of its 
Hamiltonian with reasonable computer resources. 

A few remarks are in order. 

One merit of this method is that it enables us to numerically evaluate expectation values in 
various quantum systems using relatively small portion of the whole vector space, which therefore 
enlarge the size of the systems we study. Another merit is that the method can be applied to 
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systems which have negative sign problem. Actuahy a work on the Shastry-Sutherland model^'^^ 
near the critical point as well as on the quantum spin model on a triangular lattice is in progress. 

More studies should be necessary on two issues of the method. One of them is how one can 
effectively generate the approximate state | ^a)- Although we can improve our estimation by 
increasing the number of samples ?^smpb 

it is desirable, in order to perform efficient estimations, to 
employ a good | V'a) for the following reasons. One reason is that with the better | -^a) the sample 
error will be smaller. Another reason is that with the better | -^a) the smaller L will be enough 
to obtain reliable estimates. In addition to the very empirical way we adopted to obtain | ^a) 
in this paper, a systematic way which makes use of the Lanczos method in small vector spaces 
are under investigation. There we start with a trial state | V'tri) , calculate the matrix elements of 
the Hamiltonian H, diagonalize the matrix and define | -^a) as its lowest energy eigenstate. The 
basis for the matrix elements is { | V'tri) j H \ V'tri) > H'^ \ V'tri) > • " " i \ V'tri) } with a small value of 
p, which we numerically calculate without any approximation. Another important problem is to 
find a sophisticated way to extract physical quantities from the expectation values. The results of 
several fits we tried seem to suggest there are some shortages in our simple assumptions. Both of 
the issues, which might be largely model-dependent, should be investigated in the future work. 

The on-off probability method could be applied to any problem which is mathematically related 
with the Markov process. In the Markov process, which can be described by repeated operation of a 
probability matrix Mp to an initial state v , rapid increase of states would occur. If we can replace 
(Mp)-^ V with M^jj(L)-j.Mp ■ ■ ■ M^^(i)-^Mp v the number of states we should consider becomes much 
smaller. We hope our method is helpful in many fields of science. 
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e 


(('S'{r;}))smpl 


4 


Pi 


m) 


((^))smpl 


0.3 


73.572±0.091 


21.07 


20.81 


245.24 


245.3 


0.2 


73.582±0.073 


13.76 


13.62 


366.99 


367.0 


0.1 


73.558±0.050 


6.46 


6.34 


731.97 


731.8 


0.05 


73.591±0.033 


2.84 


2.77 


1451.72 


1452.1 



Tabic I. Results on S = ^ \ci\ calculated by the on-off probability method, where Ci denotes the coefficient of state 
\i) for the exact ground state of the 4x4 Heisenberg quantum spiu system. The number of samples, Wsmpi, is 10*. 
The exact value of S is 73.57224. The values of a| and ((A'')) calculated by (8) and (9) are also shown in the table. 
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e 


((-£^E{r?}(l)))smpl 






0.20 


-11.2128±0.0215 


1.1703 


1.168 


0.15 


-ll.2198zb0.0170 


0.7345 


0.738 


0.10 


-11.2285±0.0122 


0.3606 


0.374 


0.08 


-ll.2212zb0.0095 


0.2283 


0.228 


0.05 


-ll.2230zb0.0052 


0.07338 


0.0719 


0.03 


-ll.2288zb0.0033 


0.02731 


0.0274 


0.02 


-ll.2266zb0.0021 


0.01037 


0.0106 


0.01 


-ll.2283zb0.0007 


0.00120 


0.0012 



Table II. Results on Ee{1) = {tpE\H\i'E) obtained for the 4x4 Heisenberg quantum spin system by the on-off 
probability method with nsmpi = iC, where [V'e) is the exact eigenstate of the system. The exact value is 
E = -11.2285. 



A New Approach to Stochastic State selections in Quantum Spin Systems 



13 



e 


((^E{j,} (2)))smpl 


4w(2) 




0.20 


126.200±0.344 


297.23 


296.2 


0.15 


126.006zb0.262 


170.68 


171.8 


0.10 


125.988±0.175 


77.04 


76.7 


0.08 


126.006±0.138 


47.40 


48.2 


0.05 


126.062±0.013 


14.78 


14.9 


0.03 


126.088±0.046 


5.332 


5.30 


0.02 


126.079±0.028 


1.997 


2.01 


0.01 


126.071zb0.009 


0.228 


0.23 



Table III. Results on Ee{2) = (ipElH^lipE) obtained for the 4x4 Heisenberg quantum spin system by the on-ofF 
probability method with rigmpi = 10''. The exact value is = 126.0788. 
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e = 0.1 




e = 0.01 




L 


E^{L) 










1 


-0.112285 


xlO^ 


-(0.11228 ± 0.00012) 


xlO^ 


-(0.112284 ± 0.000007) 


xlO^ 


2 


0.126079 


xlO^ 


(0.12599 ± 0.00017) 


xlO^ 


(0.126078 ± 0.000009) 


xlO^ 


3 


-0.141567 


xlO^ 


-(0.14163 ± 0.00023) 


xlO^ 


-(0.141567 ± 0.000012) 


xlO^ 


4 


0.158959 


xlO^ 


(0.15900 ± 0.00031) 


xlO^ 


(0.158952 ± 0.000016) 


xlO^ 


5 


-0.178486 


xlO^ 


-(0.17846 ± 0.00040) 


xlO^ 


-(0.178481 ± 0.000020) 


xlO^ 


6 


0.200413 


xlO^ 


(0.20023 ± 0.00051) 


xlO^ 


(0.200408 ± 0.000023) 


xlO^ 


7 


-0.225034 


xlO^ 


-(0.22498 ± 0.00067) 


xlO^ 


-(0.225029 ± 0.000028) 


xlO^ 


8 


0.252679 


xlO^ 


(0.25218 ± 0.00081) 


xlO^ 


(0.252684 ± 0.000034) 


xlO^ 


9 


-0.283720 


xlQio 


-(0.2835 ± 0.0010) 


xlOio 


-(0.283715 ± 0.000040) 


xlQio 


10 


0.318574 


xlO^^ 


(0.3182 ± 0.0015) 


xlO^i 


(0.318591 ± 0.000047) 


xlO^^ 



Table IV. Results on Ee{L) — {iPe\H'"\iPe) (L = 1,2, ■ ■ ■ , 10) obtained for the 4x4 Heisenberg quantum spin system 
by the ou-off probability method with 10"' samples. We present the data for e = 0.1 and e = 0.01 together with 
the exact values of Ee{L). 



A New Approach to Stochastic State selections in Quantum Spin Systems 



15 





e = 


0.1 


e = 


0.01 


L 


((A^,E(L))),„,pi 




((A^,^(i)))smpl 


((A^f(i)))s„.pi 


1 


12870 


731.7 


12870 


6249.41 


2 


7959.5 


564.7 


12797.06 


6242.70 


3 


6664.4 


515.6 


12796.22 


6242.22 


4 


6238.6 


498.5 


12796.22 


6242.85 


5 


6083.2 


491.9 


12796.34 


6242.79 


6 


6022.0 


488.9 


12796.25 


6243.20 


7 


5995.0 


487.9 


12796.27 


6242.38 


8 


5985.5 


487.3 


12795.97 


6243.45 


9 


5981.4 


487.2 


12796.42 


6242.71 


10 


5980.1 


487.5 


12796.21 


6243.33 



Table V. Numbers of non-zero coefficients before and after operating the random choice matrix M^^(i)j to the 
state HMf^^{L~i)^ ■ ■ ■ HM^^(i)y\tl)^)^ where H and \ipE,) denote the Hamiltonian and the exact ground state of the 
4x4 Heisenberg quantum spin system, respectively. The number of samples is 10*. 
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L 


Ea{L) 




((^A{r?}(^)))smpl 




((iVA(^)))smpl 




1 


-0.107445 


xlO^ 


-(0.10750 ± 0.00012) 


xlO^ 


5382.0 


540.97 


2 


0.119302 


xlO^ 


(0.11937 ± 0.00017) 


xlO^ 


6342.2 


466.83 


3 


-0.133162 


xlO^ 


-(0.13327 ± 0.00025) 


xlO^ 


5739.5 


445.39 


4 


0.149066 


xlO^ 


(0.14936 ± 0.00038) 


xlO^ 


5550.9 


438.58 


5 


-0.167056 


xlO^ 


-(0.16720 ± 0.00062) 


xlO^ 


5488.3 


436.49 


6 


0.187354 


xlO^ 


(0.1876 ± 0.0012) 


xlO^ 


5468.6 


435.69 


7 


-0.210193 


xlO^ 


-(0.2121 ± 0.0028) 


xlO^ 


5462.0 


435.78 


8 


0.235870 


xlO^ 


(0.2368 ± 0.0034) 


xlO^ 


5461.7 


435.47 


9 


-0.264725 


xlQi" 


-(0.2647 ± 0.0033) 


xlQi" 


5459.3 


435.24 


10 


0.297136 


xlO^^ 


(0.2958 ± 0.0049) 


xlO^^ 


5456.6 


435.58 



Table VI. Results on Ea{L) = {iljA\H^\ijjA) for the 4x4 system obtained from 10* samples with e = 0.1 and S = 0.01. 
We also present the exact values of Ea{L) in the second column for comparison. 
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L 


((^A|,/}(^)))smpl 




((iVA(L))),„,pi 




1 


-(0.235605 ± 0.000005) 


xlO^ 


1173554 


55121.7 


2 


(0.557681 lb 0.000023) 


xlO^ 


2159459.8 


72578.1 


3 


-(0.132179 ± 0.000008) 


xlO^ 


2952153.5 


88743.8 


4 


(0.313746 ± 0.000029) 


xlO^ 


3677982.2 


103182.8 


5 


-(0.74565 ± 0.00013) 


xlO^ 


4301258.0 


115698.3 



Table VII. Results on Ea{L) for the 6x6 system obtained from 5.5 x 10^ samples with e = 0.0057. The exact value 
of i;A(l) = (V-'A I H 1 -iPa) is calculated to be -23.56093. 



L 


{{EA{r,}{L)))smpl 




(«(^)))smpl 


((^a^(^)))smpl 


1 


-(0.407475 ± 0.000075) 


xlO^ 


9549922 


96800.9 


2 


(0.167062 ± 0.000054) 


xlO^ 


9572198.2 


73091.7 


3 


-(0.68731 ± 0.00047) 


xlO^ 


7155664.2 


63898.6 


4 


(0.28338 ± 0.00090) 


xlO^ 


6203627.8 


60498.0 


5 


-(0.11633 ± 0.00054) 


xlO^ 


5857166.8 


59281.9 



Table VIII. Results on Ea{L) for the 8x8 system obtained from 1.1 x 10^ samples with e = 0.01. The exact value 
of .Ba(I) is -40.74998. 
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Fig. 1. Distributions of EE{r,}^{L) (fe = 1, 2, • • • , 10'') on the 4x4 lattice for several values of L. Here the ratio 
the horizontal axis is defined as 50 times of E^{^y^{L) / E^{L) . 
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Fig. 2. Distributions of -EA{j7}k(^) = li 2, • • • , lO*) on the 4x4 lattice for several values of L. The ratio is 
times of E/^{^y^{L)/EA{L). 



